Influence of interface structure on electronic properties 
and Schottky barriers in Fe/GaAs magnetic junctions 
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The electronic and magnetic properties of Fe/GaAs(001) magnetic junctions are investigated using 
first-principles density-functional calculations. Abrupt and intermixed interfaces are considered, 
and the dependence of charge transfer, magnetization profiles, Schottky barrier heights, and spin 
polarization of densities of states on interface structure is studied. With As-termination, an abrupt 
interface with Fe is favored, while Ga-terminated GaAs favors the formation of an intermixed layer 
with Fe. The Schottky barrier heights are particularly sensitive to the abruptness of the interface. 
A significant density of states in the semiconducting gap arises from metal interface states. These 
spin-dependent interface states lead to a significant minority spin polarization of the density of 
states at the Fermi level that persists well into the semiconductor, providing a channel for the 
tunneling of minority spins through the Schottky barrier. These interface-induced gap states and 
their dependence on atomic structure at the interface are discussed in connection with potential 
spin-injection applications. 

PACS numbers: 73.20.At, 75.70.Cn, 73.30.-l-y, 72.25. Mk 
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I. INTRODUCTION 

The problem of spin injection has been a subject of in- 
tense study since the proposal of the electronic analog of 
the electro-optic modulatori Since then the entire field 
of spintronics has been developed. The main idea lies in 
the possibility of controlling the spin of charge carriers, 
thereby adding an additional degree of freedom to exist- 
ing semiconductor-based electronics. A number of novel 
electronics devices based on this idea have been proposed, 
such as reprogrammable logic devices, spin valves, spin- 
injection diodes, and devices utilizing giant magnetore- 
sistance (see, for instance, Ref. A major challenge 
in the field has been the creation of spin-polarized cur- 
rents in nonmagnetic semiconductors. One approach is 
the use of ferromagnetic contacts as spin sources. A spin 
polarization of the current is expected from the differ- 
ent conductivities resulting from the different densities of 
states for spin-up and spin-down electrons in the ferro- 
magnet. Significant progress in molecular beam epitaxy 
has allowed growth of high quality, virtually defect-free 
junctions between magnetic materials and semiconductor 
substrates, and films of ferromagnetic metals such as Fe 
or Co grown epitaxially on semiconductor structures are 
promising candidates for spin injection. 

Zhu and co-authors^ have demonstrated efficiencies of 
2% for injection of spin-polarized electrons from a metal 
into a semiconductor for a GaAs/ (In, Ga) As light emit- 
ting diode (LED) covered with Fe. Moreover, Hanbicki 
et alX have managed to achieve a spin injection efficiency 
of 30%. In the latter case, an Fe film grown on an 
AlGaAs/GaAs quantum well LED structure was used. 
In both cases it was suggested that the spin injection 
arises from tunneling of spin-polarized electrons from the 
metal into the semiconductor across the Schottky barrier 
formed at the interface. Such a tunneling process is be- 
lieved to be responsible for the spin injection since it is 



not affected by the conductivity mismatch^ between the 
metal and semiconductor that severely limits the spin- 
injection efficiency in the diffusive transport regime. A 
remarkable consequence of such a mechanism is that it 
is independent of temperature. In fact, in Refs. 3 and 
Q, nearly constant tunneling efficiencies of 2% and 30% 
were observed for a range of temperatures from 2 K to 
300 K and from 90 K to 240 K, respectively. Recently 
spin injection efficiencies of 13% have been reported at 
5 K across a Fe/GaAs(110) interface^S and 6% across 
a Fe/Ala:Gai_2;As/GaAs Schottky contact at 295 Ki 
At the moment, the room temperature record of 32% 
is held by CoFe/MgO injectors grown on p-GaAs(lOO) 
substrates 1^ 

The electronic structure of the metal/semiconductor 
interface plays an important role in spin-dependent trans- 
port properties across such junctions. For the case of 
the Fe/GaAs(001) interface, a number of first-principles 
studies have been carried out. Green's function methods 
have been used to study the electronic structure, charge 
transfer, and spin polarization in Fe/GaAs/Fe(001) tun- 
nel junctions in which the interface geometry is atomi- 
cally abrupt and ideal Aifliii It was noted that the calcu- 
lated magnetic properties at the interface are sensitive to 
the interface structure, indicating that structural relax- 
ations could be importantiiS Other studies have focused 
on the initial stages of growth of Fe on GaAs and have 
considered how structural relaxation and intermixing of 
metal and semiconductor atoms affect magnetic proper- 
ties of thin films of Fe on GaAsiiSiiii^ Still lacking is 
an understanding of how details of the atomic arrange- 
ment at the interface affect properties directly relevant to 
spin injection, such as Schottky barrier heights and the 
nature of the interface-induced states in the semiconduc- 
tor gap through which spins may tunnel. In this paper, 
we address this by considering structural models for the 
(001) interface between Fe and GaAs that allow for dif- 
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ferent degrees of intermixing and relaxation. The effects 
of interface structure on potential, charge, and magne- 
tization profiles, Schottky barrier heights, spin-polarized 
densities of states, and interface-induced gap states are 
investigated. The results are discussed in connection with 
potential spin-injection applications. 

The paper is organized in the following way. Compu- 
tational details are given in Section II. In Section III we 
describe and discuss the structures considered and the 
resulting structural, electronic, and magnetic properties. 
Conclusions are presented in Section IV. 

II. COMPUTATIONAL DETAILS 

The bulk lattice constants of bcc Fe (2.866 A) and zinc- 
blende GaAs (5.654 A) differ by almost exactly a factor of 
two. Therefore, an interface consisting of a (001) Fe slab 
placed on the (001) surface of GaAs has a lattice mis- 
match of 1.36%. The close lattice match helps keep the 
concentration of defects at such interfaces relatively low. 
In this work, interfaces were modeled using supercells 
consisting of nine layers of Fe and nine layers of semicon- 
ductor atoms (five As and four Ga, or vice versa, depend- 
ing on the semiconductor surface termination.) This is 
sufficient to ensure that the adjacent interfaces do not in- 
teract with each other, as confirmed by both total energy 
calculations and macroscopic averages of the electrostatic 
potential and charge. 

The calculations were performed with the VASP 
program^ an ab initio density functional code that uses 
the planewave pseudopotential method. The generalized 
gradient approximation was used to treat the exchange- 
correlation part of the electron-electron interactioniiS*ii 
All calculations employed ultrasoft pseudopotentialsi^ 
and planewave basis sets with a kinetic energy cut off of 
370 eV. Monkhorst-Pack mesheaiS of 8 x 8 x 4 k-points 
were used to sample the Brillouin zone, and Gaussian 
smearing of electronic states {a = 0.2 eV) was used to 
achieve faster convergence of Brillouin zone sums with 
respect to the number of k-points. 

III. RESULTS 

A. Structure 

The Fe/GaAs interface structures considered in this 
work are shown in Fig. ^ While the bare GaAs(OOl) sur- 
face reconstructs to form dimers on the surface, recent 
calculations suggest that these dimers become unstable 
upon adsorption of Fe^ Therefore, in this work, we con- 
sider 1x1 interfaces only. We assume the Fe(OOl) and 
GaAs(OOl) slabs are aligned so that in the first complete 
metal layer, sites that would have been occupied by semi- 
conductor atoms in the absence of the interface are now 
occupied by metal atoms. An abrupt interface like this 
with no intermixing of metal and semiconductor atoms 




FIG. 1: Three models for the structure of the 1x1 inter- 
face between As-terminated GaAs(OOl) and Fe(OOl) surfaces. 
Model A is an abrupt interface, model B has one layer con- 
taining both metal and semiconductor atoms, and model C 
has two intermixed layers. 

will be referred to as model A (as in Ref. I12i() . In model 
B, the bcc metal structure is partially continued into the 
first semiconductor layer so that metal atoms occupy sites 
that are normally empty in the first layer of the semicon- 
ductor. In model C, metal atoms occupy interstitial sites 
in the second semiconductor layer as well. Both As- and 
Ga-terminated interfaces can be grown experimentally, 
and therefore, both are considered here. We find that for 
the supercells of the size considered in this work, model 
C is never energetically favorable. This agrees with Ref. 
[T^ where it was found that model C has lower energy 
only for coverages of Fe not exceeding two monolayers. 
Therefore we focus on models A and B. 

The comparison of total energies of models A and B is 
not physically meaningful since the interfaces have differ- 
ent numbers of atoms. The formation energies of these 
models are, on the other hand, physically comparable 
and given by 

Eform=Et~^N,H,. (1) 
i 

Here, Et is the total energy of the supercell, Ni is the 
number of atoms of the type i in the cell, and /ii is the 
chemical potential of the i-th atom. Thus, all energies 
from the model B calculations were adjusted by twice the 
value of the chemical potential for Fe, since each super- 
cell contains two equivalent interfaces. Test calculations 
of the formation energy using supercells with one A in- 
terface and one B interface confirm that the supercells 
are large enough to ensure that the interaction between 
interfaces in negligible. 

Since electronic and magnetic properties may differ sig- 
nificantly as a result of small changes in structural pa- 
rameters, we have relaxed all the interface models con- 
sidered. In our calculations the interfaces were relaxed 
with respect to two parameters: the distance between the 
adjacent Fe and As (or Ga) layers, c?i, and the distance 
between the first two layers in the GaAs slab, c?2 (see Fig. 
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TABLE I: Interlayer relaxations, in units of the GaAs lat- 
tice constant oq. Here, Adi is the difference between the 
relaxed and ideal (0.25ao) separation of adjacent Fe and As 
(Ga) planes at the interface, and Ad2 is the change in the 
distance between the first two planes in the GaAs slab. 





As-terminated 


Ga-terminated 




Adi Ad2 


Adi Ad2 


Model A 


0.0 0.025 


0.022 0.017 


Model B 


0.017 0.068 


0.018 0.058 



The in-plane lattice constant was fixed at the bulk 
GaAs value, ag = 5.654 A. The results of the relaxation 
are presented in Table|l] For As-terminated model A, the 
distance between As and Fe planes remains at the ideal 
value of 0.25ao, while the separation between the first 
plane of As and the adjacent plane of Ga is stretched 
by 10%. With Ga-termination, both di and ^2 increase 
by similar amounts in model A. In model B, which has 
the intermixed layer, the first two planes of the semi- 
conductor slab are repelled considerably farther apart, 
regardless of the termination. It was previously pointed 
out that lowering the concentration of atoms in the in- 
terface region is energetically favorable in model B, since 
the electrons from the extra Fe atoms in the intermixed 
layer fill antibonding orbitals and weaken the interface 
bondingpi^ 

The results of the formation energy calculations for the 
ideal and relaxed interfaces are summarized in Table ITTl 
For both ideal and relaxed geometries, model A is ener- 
getically favored for As-terminated interfaces, and model 
B is favored for Ga-terminated interfaces. This results 
from an interplay between optimization of the coordi- 
nation of the metal and semiconductor atoms and the 
relative strengths of the metal-cation and metal-anion 
bonds. Because of stronger pd hybridization, the Fe-As 
bond is more stable than the Fe-Ga bondr^^ For model 
A, the interface Fe site is six-fold coordinated, with four 
Fe neighbors in the second metal layer and two semicon- 
ductor neighbors in the first semiconductor layer. While 
putting Fe atoms in interstitial sites in the first semi- 



TABLE II: Formation energy differences, in eV per 1x1 in- 
terface unit cell. With As-termination, the abrupt interface 
of model A is favored, while with Ga-termination, the inter- 
mixed interface of model B is preferred. 





As-terminated 


Ga-terminated 




Ideal - Relaxed 


Ideal - Relaxed 


Model A 


0.046 


0.059 


Model B 


0.282 


0.291 




Model A - Model B 


Model A - Model B 


Ideal 


-0.400 


0.318 


Relaxed 


-0.164 


0.551 



conductor layer fully coordinates the interface Fe sites 
in model B, it also significantly weakens the bonding 
between the first two semiconductor layers, which be- 
come overcoordinated. At the As-terminated interface, 
the strong Fe-As bonds compensate for the undercoordi- 
nation of the Fe sites, making model A favorable, while 
at the Ga-terminated interface, the weak interaction be- 
tween Fe and Ga favors full coordination of Fe interface 
sites, as in model B. 

B. Schottky Barriers and Electronic Structure 

Figure El shows the calculated site-projected densities 
of states (DOS) for atoms in different layers in the relaxed 
As-terminated model A. The general features are similar 
for all the structural models considered. The DOS in the 
most bulk-like layers of the supercell closely resembles the 
DOS of bulk Fe or GaAs. Near the interface, a peak de- 
velops in the DOS at the Fermi level. The peak is largest 
in the first Fe layer and decreases into both the GaAs 
and Fe slabs. In the first few GaAs layers, states spread 
throughout the entire band gap of bulk GaAs, though 
the gap is practically recovered in the layer farthest from 
the interface. BardeenSi estimated that surface states 
with density ^ 10^'^ eV~^ cm~^ would effectively pin 
the Fermi level close to the charge neutrality level in the 
semiconductor gap. Our calculations for the relaxed As- 
terminated model A, for example, yield values of DOS 
at the Fermi level of ~ 2 x 10^^ eV~^ cm^^ in the first 
Fe layer, and ^ 1 x 10^'* eV~^ cm~^ in the first two lay- 
ers of the GaAs slab. These may be sufficient to pin the 
Fermi level, which consequently affects the Schottky bar- 
rier height (SBH). We will return to examine the nature 
of these mid-gap states, their spin polarization, and their 
role in spin tunneling in section D. 

As mentioned, Schottky barriers may be a necessary 
mechanism for overcoming the conductivity mismatch for 
the injection of spin from a ferromagnetic metal into a 
semiconductor. Schottky barrier heights can be evalu- 
ated from first principles using a macroscopic average 
method combined with the supercell approachi2^ Two 
necessary conditions are: 1) the supercell contains two 
equivalent interfaces, which eliminates any electric fields 
that might be present due to unbalanced charges, and 2) 
the supercell is sufficiently large so the bulk charge and 
potential properties are recovered in the most bulk-like 
layers of the supercell. In such a way an isolated interface 
is accurately modeled using a supercell. 

The calculation of the p-type Schottky barrier height, 
(j)p, is split into two parts: 

(l)p = AK„ -I- AV. (2) 

The band structure term, 

AE,^Ef-E,, (3) 

is the difference between the Fermi level Ep in the metal 
and the valence band edge in the semiconductor Ey, 
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FIG. 2: Site-projected densities of states for atoms in different 
atomic layers of the relaxed GaAs/Fe(001) supercell with As- 
terminated model A interfaces. The bottom panel represents 
the semiconductor layer farthest from the interface (As #5), 
and the top panel represents the most bulk-like layer of Fe 
(Fe #5). The vertical line at zero energy indicates the Fermi 
level. In each panel, the vertical scale runs from to 1.8 
states / eV / atom. 



where each is measured with respect to the average elec- 
trostatic potential in the corresponding bulk material. 
The band structure term is calculated from separate bulk 
calculations for the two constituents of the Schottky con- 
tact. This term implicitly includes all quantum mechan- 
ical effects as well as the exchange-correlation part of the 
potential. 

The other contribution to (j)p is the potential line-up 
across the interface AT^. This potential line-up is related 
to the dipole moment of the charge profilers depends 
on the structure of the interface, and therefore cannot 
be calculated simply from two bulk calculations. It is 
the difference between the macroscopic averages of the 
electrostatic potential in two bulk-like regions of the su- 
percell. As an example. Fig. |3| shows the planar average 
over the x — y plane and the macroscopic average of the 
electrostatic potential computed for the superlattice with 
As-terminated relaxed model A interfaces. The macro- 
scopic average lacks the bulk-like oscillations present in 
the planar average, and thereby allows one to extract the 
desired potential line-up Ay. The macroscopic average, 
however, does not provide reliable information about po- 
tential behavior in the close vicinity of the interface since 
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FIG. 3: Planar and macroscopic averages of the electrostatic 
potential calculated for the relaxed As-terminated model A 
supercell. The potential is plotted as a function of z, the 
position of atomic layers along the axis of the supercell. The 
potential line-up AV needed for the Schottky barrier height 
calculation is given by the difference between the two plateaus 
in the macroscopically averaged potential. 



the resolution of the method is limited to the lattice pe- 
riod, the GaAs lattice constant ao in our case. There is 
a notable difference between the planar average and the 
macroscopic average in the vicinity of the interface. 

We have calculated the p-type SBH for the ideal and re- 
laxed structures described in the previous section. Since 
it is electrons that are injected from Fe into GaAs in spin 
injection experiments, the n-type SBH is of more inter- 
est. To obtain the n-type SBH from our calculations, 
we subtract the calculated p-type SBH from the exper- 
imental GaAs band gap of 1.52 eV.^^ The results are 
presented in Table Hill The calculated SBH is sensitive 
to structural changes, and the differences between the 
SBH values in Table IIIII show the magnitude of changes 
one can expect for different interface structures. Inter- 
mixing of metal and semiconductor atoms at the inter- 
face decreases the n-type Schottky barrier height signif- 
icantly (by about 0.2-0.3 eV), while the effect of inter- 
layer relaxation is generally weaker. This level of sen- 
sitivity of SBH values to interface structure is consis- 



TABLE III: Calculated n-type Schottky barrier heights (in 
eV) for ideal and relaxed geometries. Interfaces with inter- 
mixing (model B) have lower barrier heights than abrupt in- 
terfaces (model A). 





As-terminated 


Ga-terminated 




Ideal Relaxed 


Ideal Relaxed 


Model A 


0.87 0.82 


1.01 1.08 


Model B 


0.58 0.64 


0.69 0.89 



tent with the observed dependence of measured SBHs 
on growth conditions. For Fe/GaAs interfaces fabri- 
cated by metal evaporation in ultrahigh vacuum, SBHs 
in the range of 0.72-0.75 eV have been reported.^'* More 
recent experiments on As-terminated atomically clean 
(001) interfaces grown by molecular beam epitaxy have 
yielded barriers around 0.90-0.92 eVi2SiS& In the latter 
work, the samples had sufhciently low interface defect 
concentrations that Fermi-level pinning by antisite de- 
fects was suppressed Our calculated results for the 
energetically favored geometries, i.e. As-terminated re- 
laxed model A and Ga-terminated relaxed model B, are 
close to the SBH values obtained for these clean samples. 
Overall, the present results compare more favorably to 
experiments than earlier density-functional calculations 
of the p-type SBH at ideal Fe/GaAs (110) interfacesiSi 



C. Charge Distribution and Magnetic Moments 

To investigate the redistribution of charge that gives 
rise to the Schottky barrier, we have computed the charge 
within Wigner-Seitz spheres centered on atomic sitesi^ 
We have also compared the planar-averaged charge in 
the interface region to the planar-averaged charge in 
the bulk-like regions and to the planar-averaged super- 
position of atomic charges. The qualitative features 
of the charge redistribution deduced from the sphere 
charges are consistent with those suggested by the planar- 
averaged charge. Figure ^ shows the difference between 
the charge inside the spheres in the supercell and in the 
bulk. At the abrupt interface (model A), there is an ev- 
ident transfer of charge from the interface Fe layer into 
the semiconducting slab, regardless of the semiconductor 
termination. In Ref. it was found that for unrelaxed 
abrupt interfaces, the interface Fc loses more charge to 
Ga neighbors than to As neighbors. We find that this 
difference is enhanced when the interlayer distances are 
allowed to relax. Interface Fe atoms also lose significant 
charge to the interstitial region where bond formation 
takes place. The redistribution of charge is limited to re- 
gions close to the interface. By the third semiconductor 
layer, the charge within the spheres has nearly recovered 
to bulk values. 

In model B, Fe sites in the first full Fe layer are fully 
coordinated, so the local charge around these sites is 
much closer to the bulk value than in model A. The 
charge redistribution takes place primarily within the in- 
termixed layer, where the Fe sites are undercoordinated 
and the semiconductor sites are overcoordinated. While 
the amount of local charge lost by Fe sites in the in- 
termixed layer is about the same for both terminations 
(~ 0.2e), with Ga-termination, the Ga sites in the in- 
termixed layer gain electrons while with As-termination, 
the As interface sites lose electrons, indicating a transfer 
of charge into interstitial regions in the Fe-As intermixed 
layer. 

Fig. [SI shows the local magnetic moments across the 
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FIG. 4: Difference between the charge inside the atomic 
spheres for the supercell and bulk GaAs and Fe, plotted as a 
function of distance z along the direction normal to the inter- 
face. The upper panels show results for relaxed As-terminated 
models A and B, and the lower panels show results for relaxed 
Ga-terminated models A and B. In the intermixed layer in 
model B, the higher point represents the charge difference for 
the As or Ga atom and the lower point represents that for the 
Fe atom. 




FIG. 5: Local magnetic moments plotted as a function of 
distance z along the direction normal to the interface. The 
upper panels show results for relaxed As-terminated models 
A and B, and the lower panels show results for relaxed Ga- 
terminated models A and B. In the intermixed layer in model 
B, the higher data point represents the Fe moment while the 
lower data point represents the Ga or As moment. 



Fe/GaAs interface for all four structurally relaxed mod- 
els. The magnetic moments were calculated from the 
integrated spin-polarized DOS within the Wigner-Seitz 
spheres. There is a small enhancement of the magnetic 
moment near the Fe interface in comparison to the bulk- 
like moment of '^2.3/iB in the central Fe layer. This is 
similar what was found in Ref.^Jl- Since the Fe mo- 
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FIG. 6: Local majority-spin (circles) and minority-spin (di- 
amonds) densities of states at the Fermi level in different 
atomic layers. The top panels show results for relaxed As- 
terminated models A and B, and the bottom panels show 
results for relaxed Ga-terminated models A and B. The ver- 
tical scale, in units of states/eV/atom, is kept the same in all 
panels to facilitate comparisons. 



ment at the interface is sensitive to the Fe-As(Ga) bond 
length and can be quenched by reducing the bond length 
by a few percent ji^ the relaxation of the interlayer dis- 
tances near the interface is important. In model A the 
enhancement of the Fe moment at the interface ranges 
from O.l/is for As termination to O.SfiB for Ga termina- 
tion. Model B has a larger enhancement of the magnetic 
moment at the Fe sites located in the intermixed layer, 
with spin moments ~ OAfiB larger than at the bulk-like 
sites. These results, which indicate that Fe is ferromag- 
netic at the GaAs interface, are consistent with recent 
experiments on the Fe/GaAs (100)-4x6 interface^^ where 
bulk-like spin magnetic moments were observed using x- 
ray magnetic circular dichroism. In our calculations, all 
four structural models also have small induced opposite 
magnetic moments in the first semiconductor layer with 
the approximate values of 0.04/iB and 0.02/is for As- 
terminated models A and B, respectively, and O.OSfis 
and 0.09/xb for Ga-terminated models A and B, respec- 
tively. 



D. Spin-Polarization and Interface States 

Figure|n|shows the local majority (spin-up) and minor- 
ity (spin-down) densities of states calculated at the Fermi 
energy for the relaxed As- and Ga-terminated models A 
and B. Each point on the curves corresponds to an atomic 
layer. In all four cases the majority-spin DOS is largest 
at the center of the Fe slab and decays monotonically 
towards and across the interface into the GaAs slab. At 
the same time the minority-spin DOS exhibits a sharp 
peak in the vicinity of the interface on the Fe side and 



decreases into the center of the Fe slab. In the GaAs 
slab both majority and minority DOS decay exponen- 
tially, but the dominance of the minority states is pre- 
served throughout the GaAs slab. Hence, close to the 
interface on the Fe side, the spin polarization reverses 
sign and peaks because of the large difference between 
spin-up and spin-down density of states. While model B 
exhibits such a reversal with relatively modest differences 
between spin-up and spin-down DOS, model A shows sig- 
nificantly larger values of spin polarization. 

Similar behavior of the spin-dependent DOS at the 
Fermi level has been observed both theoretically and ex- 
perimentally for the free Fe surfacei^2i2i A peak in the 
scanning tunneling spectra of the Fe(OOl) surface was at- 
tributed to a minority-spin surface band located about 
0.3 eV above the Fermi level at the F point in the two- 
dimensional Brillouin zone. We find analogous behavior 
at the Fe/GaAs interface, with the increase in minority 
spin DOS at the Fermi level at the interface attibutable 
to d states localized on Fe interface sites. The valence 
charge density for one of these interface states at the As- 
terminated abrupt interface is shown in Fig. This 
d322_r2-derived Shockley-like state, located about 0.4 eV 
above the Fermi level at the F point, produces high den- 
sities of states within the GaAs band gap and decays 
evanescently into the GaAs slab. At the intermixed in- 
terfaces of model B, analogous states localized on the 
Fe sites in the intermixed layer are found near the Fermi 
level, but because of the stronger infiuence of the reduced 
symmetry of GaAs compared to Fe, some of these inter- 
face states are more clearly a mixture of d^z^^r^ and d^y 
character, with lobes pointing along the z and x — ±y 
directions. (The x and y directions are along the cubic 
axes of both the bcc Fe and zinc-blende GaAs lattices.) 

These interface states likely play a role in Fermi-level 
pinning, which experimentally manifests as an insensitiv- 
ity of the SBH to the metal workfunction. Our calcula- 
tions suggest that at the defect-free interface, the Fermi 
level is pinned by Fe minority-spin interface states, sup- 
porting the metal-induced gap states (MIGS) model22i2^ 
rather than the semiconductor surface state modeliSi 
This also agrees with recent experiments on the pressure 
dependence of metal/GaAs Schottky barrier heightS)2& 
which support the MIGS model for atomically clean 
Fe/GaAs interfaces. A notable difference between the 
interface states we find and MIGS suggested in the orig- 
inal works^SiS, is that the interface states in our case are 
localized at the interface and are not derived from the 
metal bulk states. 

To quantify the decay of the spin polarized DOS into 
the GaAs, we have performed additional large supercell 
calculations using the relaxed interface geometries, in- 
creasing the number of atomic layers in the GaAs slab 
to 17. A large number of k-points was used in the DOS 
calculations (24x24x4), and a Gaussian broadening of 
states with ct = 0.1 eV was used. We fit the DOS at 
the Fermi level in the GaAs slab to the functional form 
g~2K2^ where k is a decay constant, and z is the distance 
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FIG. 7: Valence charge density of a minority-spin interface 
state at the F-point in the relaxed As-terminated model A 
supercell. The top half of the figure contains Fe layers and 
the bottom half contains GaAs layers. This state lies about 
0.4 eV above the Fermi level. 
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FIG. 8: The site-projected GaAs DOS at the Fermi level as a 
function of the distance from the Fe slab. The majority-spin 
DOS are plotted with circles and the minority-spin DOS are 
plotted with triangles. The upper panels correspond to As- 
terminated models A and B, and the lower panels correspond 
to Ga-terminated models A and B. The DOS are in units of 
states / eV / atom. 



from the interface. The decay constants k-^ j for states 
with different spin, and corresponding decay lengths 
are listed in Table IIVI Semi-logarithmic plots of the 
majority- and minority-spin DOS at the Fermi level in 
the GaAs slab are shown in Fig. |S1 along with the fit- 
ted exponentials (dashed lines). Each point corresponds 



to an atomic layer, and data corresponding to the in- 
termixed layer in model B was omitted since this layer 
contains both Fe and semiconductor atoms. At the Fermi 
level the minority-spin states continue to dominate the 
majority-spin states throughout the GaAs slab. There is 
a significant proximity effect in the sense that the ratio 
of the DOS corresponding to the two spin states does 
not change significantly. In other words, the spin polar- 
ization of the states at the Fermi level remains virtually 
constant into the bulk of GaAs, even as the magnitude 
of the DOS decays exponentially. 



TABLE IV: Decay constants (in A ^) for majority and 
minority states and corresponding decay lengths (in A). 





As-terminated 


Ga-terminated 












Model A 


0.17 


0.16 


0.15 


0.16 


Model B 


0.16 


0.14 


0.15 


0.18 




h 




't 




Model A 


3.0 


3.1 


3.3 


3.1 


Model B 


3.2 


3.5 


3.4 


3.0 



The calculated decay constants k^^i listed in Table Hvl 
show that there is not a significant difference between the 
decay of spin-up and spin-down states for a given struc- 
ture geometry. Such behavior is expected (assuming that 
the effective masses of spin-up and spin-down electrons 
do not differ significantly) since states of the same en- 
ergy should have the same decay constants unless the 
potential barriers for the two particles differ. The simi- 
larities of decay constants show that the Schottky barrier 
heights for both spin-up and spin-down electrons are the 
same, which is reasonable since the GaAs slab is essen- 
tially non-magnetic. This finding, however, is in contrast 
with results of calculations of the electronic properties of 
C0/AI2O3/C0 magnetic tunnel junctions, where differ- 
ent decay constants for the spin-polarized states at the 
Fermi energy in AI2O3, and therefore different potential 
barrier heights, were obtained for electrons of opposite 
spins. In terms of structure dependence, we find that the 
interfaces with intermixing of metal and semiconductor 
atoms tend to have longer decay lengths than abrupt in- 
terfaces, which is consistent with trends in the calculated 
Schottky barrier heights. 

While detailed transport calculations, using, for ex- 
ample, the Landauer or Kubo formalisms, are needed 
for quantitative predictions regarding spin-polarized cur- 
rent, symmetry considerations allow us to make some 
qualitative arguments about the spin-injection process 
through the interface. Electrons tunnel from bulk states 
near the Fermi level in Fe through evanescently decaying 
gap states in the Schottky barrier region into states near 
the conduction band minimum in bulk GaAs (Fig. O. 
The states participating in this transport process must 
be compatible by symmetry. In the approximation that 
transport across the interface is dominated by carriers 
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FIG. 9; Schematic flat-band diagram showing that electrons 
from bulk states near the Fermi level in Fe tunnel through 
decaying gap states in the Schottky barrier into states near 
the conduction band minimum in GaAs. (Note that the width 
of the barrier region is typically much larger than the decay 
length of the gap states.) In both majority and minority spin 
channels, there are bulk-like states in each region that are 
compatible by symmetry. Only in the minority spin channel 
are there interface-derived gap states of the right symmetry 
near the Fermi level. 



with wave vector perpendicular to the interface, we fo- 
cus on states with k|| w (i.e., along the A direction in 
the cubic Brillioun zone of both GaAs and Fe). The 1x1 
interface considered here has C2v symmetry, and so does 
the lowest conduction band of GaAs (indexed as Af""^^). 
In the Schottky barrier region, the most important bulk 
states in the complex band structure are those with the 
longest decay length. In GaAs, these states have Af"-"^'' 
symmetry as welL^'^ Therefore, if we consider only bulk- 
like states, carriers must originate from states in the bulk 
Fe bandstructure that are compatible by symmetry with 
^GaAs s|;a,tes. In the majority-spin band structure of Fe, 
the (i322_r2-derived band, indexed as Af^, crosses the 
Fermi level and satisfies the symmetry requirement. In 
the minority-spin band structure, the (ia:j,-derived band, 
indexed as Af^*^, crosses the Fermi level and satisfies the 
symmetry requirement. Since the Af"^ band has some 
s character as well as d^^2_j.2 character, it is expected 
that the coupling of states across the interface will be 
stronger in the majority channel. Thus, considering only 
bulk-like states, the majority spin current would be ex- 
pected to dominate. In fact, transport calculations based 
on the Landauer formalism support this conclusion. '^^ 

However, in addition to bulk-like states, interface 
states may play a role in the tunneling process. The 
interface states in the vicinity of the Fermi level all have 
minority spin. Arising from d^^^^j-i and dxy orbitals on 
interface Fe sites, and resonant with the bulk Fe A^,'^ 
minority band that crosses Ep, these states provide ad- 
ditional symmetry-compatible gap states for tunneling of 
minority spins through the barrier. Hence the presence 



of these interface states close to the Fermi level could 
reduce the spin polarization of the tunneling current, or 
even reverse its sign. In Tabled we list the energy of 
the symmetry-compatible interface states within 0.5 eV 
of the Fermi level at the F point for different structural 
models. These states typically have dispersions of a few 
tenths of eV across the Brillouin zone. Since the prox- 
imity of these interface states to the Fermi level depends 
on the atomic structure at the interface, we expect the 
spin polarization of the injected current to be sensitive 
to interface structure as well. Measurements of the sign 
of the circular polarization of electroluminescence from 
Fe/GaAs spin LEDs indicate injection of minority spins^ 
opposite to what is predicted from transport calculations 
for ideal junctions.-^ The present results indicate that 
details of the interface structure, such as the degree of 
intermixing and relaxation, likely contribute to this dis- 
crepancy. 



TABLE V: Energy of symmetry-compatible minority-spin in- 
terface states within 0.5 eV of the Fermi level. Energies are 
in eV and measured relative to the Fermi level. 





As-terminated 


Ga-terminated 




Ideal Relaxed 


Ideal Relaxed 


Model A 
Model B 


0.15, 0.16, 0.37 0.31, 0.39, 0.39 
0.34, 0.41 0.16, 0.48 


N/A 0.43 
-0.02 0.12, 0.28 



IV. CONCLUSIONS 

We have investigated the electronic and magnetic prop- 
erties of Fe/GaAs (001) interfaces using density func- 
tional calculations. Two structural models were com- 
pared: one with an abrupt interface and one with in- 
termixing of metal and semiconductor atoms. Both As- 
and Ga-terminations were considered, and interlayer sep- 
arations were relaxed. Due to differences in Fe-As and 
Fe-Ga bonding, the As-terminated structure favors the 
abrupt interface while the Ga-terminated structure fa- 
vors intermixing at the interface. Magnetization profiles 
show bulk-like magnetic moments at the interfaces. In 
all cases, charge is transferred from Fe to GaAs, creating 
Schottky barriers that vary in height depending on de- 
tails of the interface structure. In general, the SBHs are 
less sensitive to interlayer relaxations than to the nature 
of the interface, with intermixing of atoms at the inter- 
face leading to smaller n-type barrier heights. In all the 
structural models considered, the minority-spin Fe inter- 
face state induces states of Af"'^'' symmetry within the 
semiconductor gap. These states lead to a reversal of the 
sign of the spin polarization of the density of states near 
the interface (compared to the Fe bulk), and this spin 
polarization of the density of states at the Fermi level 
persists well into the semiconductor. These interface- 
induced gap states (which also are likely to play a major 
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part in pinning of the Fermi level) provide an additional 
channel for tunneling of minority spins. The proximity 
of these interface states to the Fermi level, which affects 
the magnitude and possibly the sign of the spin polar- 
ization of the tunneling current, varies significantly with 
interface structure. 
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